Properties of the quark-antiquark-gluon Fock state in the rj c meson 
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In this paper the twist-3 distribution amplitude of the quark-antiquark-gluon Fock state in the 
n c meson is studied. To calculate the moments of this distribution amplitude QCD sum rules is 
applied. Using the results of the calculation the model of this distribution amplitude is built. In 
addition NRQCD matrix elements which determine the properties of the quark-antiquark-gluon 
Fock state are determined. In particular, the probability amplitude to find the quark-antiquark pair 
in the color octet state, the mean gluon energy and the fraction of momentum carried by gluon, the 
relative velocity of the color-octet quark-antiquark-pair are calculated. 

PACS numbers: 12.38.-t, 12.38.Bx, 

I. INTRODUCTION. 



Theoretical description of hard exclusive processes is based on the factorization theorem [l|, Q ■ Within this theorem 
the amplitude of hard exclusive process can be separated into two parts. The first part is partons production at very 
Q\ ' small distances, which can be treated within perturbative QCD. The second part is hardronization of the partons at 
CS| ■ larger distances. This part contains information about nonperturbative dynamic of the strong interactions and it can 
be parameterized by process independent distribution amplitudes (DA), which can be considered as wave functions 
pH . of hadrons at light-like separation between the partons. 

Recently, two-particle charmonia DAs have become an object of intensive study 04I3|- Research in this direction 
allowed one to build models for two-particle charmonia DAs and to carry out the study of some interesting exclusive 
processes with charmonia production (see, for instance, [l5l - [l7j ). 

The first aim of this paper is to study the properties of the twist-3 quark-antiquark-gluon DA of the r\ c meson within 
i— QCD sum rules and to build the model of this function which can be used in the calculation of different production 
processes. There are a lot of papers (see, for instance, 0, 0, [Hj]) devoted to the study of the quark-antiquark-gluon 
Fock state DAs of light pseudoscalar mesons. However, until now there are no papers devoted to the study of similar 
■ DAs for charmonia. 

\q \ Another description of charmonia production processes, which is called nonrelativistic QCD(NRQCD), is based on 
f^) ■ the expansion of matrix elements in powers of the relative velocity of quark-antiquark pair inside charmonia mesons 
^sO \ [20| . NRQCD matrix elements play crucial role in this approach. These matrix elements parameterize charmonia 
structure. For instance, there is the matrix element which determines the probability amplitude to find quark- 
antiquark pair in the color-octet state in the rj c meson. Another example is the matrix element which determines the 
gluon energy in the r\ c meson. The second aim of this paper is to determine some of these matrix elements from the 
11 quark-antiquark-gluon DA of the r\ c meson. 

This paper is organized as follows. In the next section the definition of the twist-3 quark-antiquark-gluon DA 
will be given. The parameters of this DA will be calculated in the third section. The forth section is devoted to 
the calculation of the NRQCD matrix elements which parameterize the quark-antiquark-gluon Fock state in the rj c 
meson. In section V the model of the DA under consideration will be proposed. In the last section the results of this 
paper will be summarized. 

II. DEFINITION. 

The twist-3 distribution amplitude (DA) of the quark-antiquark-gluon Fock state in the r\ c meson can be defined 
as follows [3 Q1 

(Vc(p)\q{z)[z 1 vz}a^ 5 gG a j 3 {vz)[vz 1 -z}q(-z)\0)Q = ifz v {p a Pt,9ufs ~ P<xP»9»p ~ PpP»9va + PpPvQ'^x) x W 

x J dx 1 dx 2 dx 3 S{x 1 +X2+x 3 -l)e i ^ Xl - X2+VX3) V Q (x 1 ,X2,x 3 ), 
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where Vq(xx, X2, x 3 ) is the twist-3 DA of the quark-antiquark-gluon Fock state, xi, X2, x 3 are the fractions of momen- 
tum of the rj c meson carried by the quark, antiquark and gluon correspondingly, 

9^ - 9tiu — , 

[zx ,z 2 ]=P exp (ig J dz^A^j . 

Below it will be assumed that the function Vq(xi, X2, x 3 ) is normalized as follows 

dxidx2dxz8(xi + x 2 + x 3 - 1)Vq(xi, x 2 , £3) = 1 (2) 



The DA Vq(xi, X2, X3) is defined at the renormalization scale Q. The renormalization properties of three-particles 
DAs are rather complicated and they will not be discussed in this paper. For this reason below the subscript Q will 
be omitted. It is worth to note that if the energy scale tends to infinity Q — > 00, the DA Vq(x%, X2, X3) tends to it's 
asymptotic form 

V as (xi, x 2 , x 3 ) = 360x1x2X3. (3) 

Because of nonrelativistic nature of heavy quarkonia, one can expect that real DA is very far from its asymptotic 
form. 

Another very useful property of the DA is that the function V(xi,X2,X3) is symmetric under the replacement 
xi O X2: 

V(x 1 ,x 2 ,x 3 ) = V(x 2 ,x 1 ,x 3 ). (4) 
Below we will also need the relation between the moments of the DA V(x%, X2, X3) and QCD matrix elements 

2i/ 3 ^(^)" 1+n2+n3+2 K 1 < 2 ^3 3 } = <r?c(p)|{^Sr}a M75 { 5 (i^)"«G^4{(^r 9 }|0> Z V. (5) 



The properties of the DA V(xi, X2, x 3 ) can be parameterized by the moments (x^x^Xg 3 ). In this paper the following 

4) 



moments will be calculated (X3), (x 3 ), ((xi — X2) 2 ). The other the first and the second moments can be expressed 



through these ones. 



III. THE MOMENTS OF THE DISTRIBUTION AMPLITUDE V{xi,x 2 ,x 3 ) FROM QCD SUM RULES. 

QCD sum rules [2ll [22| proved to be very effective and universal tool in the determination of different matrix 
elements. In particular,QCD sum rules can be applied to study the properties of DAs. There are a lot of papers 
( see, for instance, @, [3 ) where QCD sum rules was applied to study the parameters of the twist-3 DAs 
of light pseudoscalar mesons. However, there are no papers devoted to the study of heavy quarkonia DAs of the 
quark-antiquark-gluon Fock state. 

In this section QCD sum rules will be applied to calculate the moments of the function V(xi,X2,x 3 ). To do this 
let us consider the correlator 

n ni , n2 , n3 (q,z) = i J d 4 xe^(0|TJ+, l2j „ 3 (x,z)J ,o,o(0,z)|0) = (^)™ 1+ ™ 2 +" 3+4 n(g 2 ), (6) 

where the current J ni ,n 2 .n 3 (x, z) is defined as follows 

J ni ,n 2 ,n s {x,z) = {g(izS)" 1 } ^ 7 AlA 75{ 5 (iz5)' l3 G All ,}{(^z5)" 2 (Z }z A z^ (7) 
Applying standard procedure which will not be described here one gets QCD sum rules for correlator (J6j) 
lf2 (x?x?z?) _1 [ S0 . Imn pcrt ( S ) 1 / d \ m 

where Q 2 = —q 2 , n pcrt (s), n nport (s) are the perturbative and nonperturbative contributions to correlator ©. The 
nonpertubative part of the correlator parameterizes the contribution of the vacuum condensates, which at the leading 
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order approximation used in this paper is given by the gluon vacuum condensate (a s G 2 ). The expressions for the 
ImIIp Crt (s), n npert (s) are rather lengthy and very complicated. For this reason they will not be shown here. 

In the original paper [22[ the method QCD sum rules was applied at Q 2 = 0. However, as was shown in paper [23} 
there is a large contribution of higher dimensional operators at Q 2 = which grows rapidly with m. To suppress this 
contribution in this paper sum rules ((SJ will be applied at Q 2 = 4m 2 . 

In the numerical analysis of QCD sum rules the values of the parameters m c and (a s G 2 ) will be taken from paper: 
d 

m c = 1.24 ±0.02 GeV, (— G 2 ) = 0.012 ± 30% GeV 4 . (9) 

The value of the threshold So is varied within the interval ^/sq € (3.7 ± 0.2) GeV. The value of the strong coupling 
constant is taken from paper [24[ 

a s {m T ) = 0.332 ±0.015, (10) 

and then it is evolved to the scale [i — m c at two loops. In the calculation the uncertainty of the results due to the 
variation of the strong coupling constant a s is not very important. For this reason it will be disregarded. 
The application of sum rules (j8]) to the calculation of the constant f 3v gives the following result 

flj = (12±3±2± 1) • 10" 6 GeV 4 , (11) 

where the second, the third and the forth errors are due to the uncertainties in the values of the parameters m c , 
(a s G 2 ) ((9]) and so correspondingly. 

Further let us proceed to the calculation of the moments (x™ 1 x 2 2 x 3 3 ) . The result of the calculation of the moment 
(x 3 ) is 

(x 3 ) = 0.22 ± 0.01 ± 0.01 ± 0.02. (12) 

The sources of errors are the same as in (fTTj) . Evidently, due to the symmetry relation (|4|) the other hrst moments 
can be related to the (23) as follows (xi) = (x 2 ) = ((1 — x 3 )/2). 

It should be noted that moment (II 2[) has simple and very important physical meaning. The (x 3 ) measures the 
fraction of momentum carried by gluon in the quark-antiquark-gluon Fock state of the rj c meson. Value (|12p tells us 
that the gluon carries ~ 20% of the total momentum. 

Another moments which also have important physical meaning are ((xi — X2) 2 ), (x 2 ). The calculation gives 

((11 - x 2 f) = 0.021 ± 0.002 ± 0.003 ± 0.001 

(xj) = 0.055 ± 0.005 ± 0.006 ± 0.004. (13) 

The sources of errors are the same as in (fTT|) . It is easy to show that all moments of the first and the second order 
can be expressed through moments (fl"2"|) . (fT3"|) . 

It should be noted here that matrix elements (|TT1) . (fT2|) . (fl3l) are scale dependent quantities. The virtuality of the 
propagators in the calculation of the ImiI por t(s), n npC rt(s) is of order of ~ m c . So, results (fTT|) . (fT2|) . (Q~3|) are defined 
at the same scale. 



IV. NRQCD MATRIX ELEMENTS AND PROPERTIES OF THE QUARK-ANTIQUARK-GLUON 

FOCK STATE IN THE Vc MESON. 

As was noted in the introduction, the DA V(xi, X2, x 3 ) contains the information about properties of the quark- 
antiquark-gluon Fock state in the rj c meson. According to formula (J5]) this properties can be parameterized by some 
rather complicated QCD operators. To make the result of this paper more transparent, one can expand these operators 
in the relative velocity of the quark-antiquark pair in the i] c meson, as was done in paper Q . Following the approach 
proposed in this paper, at the leading order approximation one gets the relations 

(^Mffi^xlO) - 3f 3v M 2 , (14) 

(»fc|^ + ffi((«Vo)$ffi)x|0) = 3/3,^(13), (15) 

(t?c|^ + o-i((-3Vo - V^OxlO) = 9f 3v M*(x 2 ) 7 (16) 

( Vc \^ + <T 3 ((-V 2 k gH 3 - gH j V k + 2V k gH J V k )) X \0) = 9f 3n M*((x 1 - x 2 ) 2 ). (17) 



4 



where ip + and x are Pauli spinor fields that create a quark and an antiquark respectively, a is the Pauli matrices, 
H is a chromomagnetic field, j, k are spatial indexes. Note that in (jT4"]) - (fTT|) summation over the repeated indexes is 
assumed. 

From formulas (|14[) - (fT7)) one sees that the quark-antiquark pair is in the S-w&ve spin-triplet color-octet state. The 
role of the color-octet state in the r/ c meson is determined by the matrix element 

(77 c |^ + 3 (ai?) X |0) = 3f 3n M 2 = (92 ± 34) • 1CT 3 GeV 4 . (18) 

Note that in the last formula the uncertainty due to the higher order terms in the relative velocity expansion have 
been taken into account. 

Strictly speaking, the role of the color-octet state in the rjc meson strongly depends on the process. However, one 
can estimate the probability amplitude to find the color-octet state as follows 

(qqG\r/ c } 0.04. (19) 

One sees that this value is very small. In paper [2"j the ratio (qqG\r) c ) was estimated as ~ v 7 / 2 ~ 0.06 — 0.09 for 
v 2 ~ 0.20 — 0.25, what is in reasonable agreement with result (fH)|) . 

Now one can estimate the average fraction of momentum carried by gluon in the r\ c meson. To do this one can 
apply the standard quantum mechanical formula 

{Xg)r, = ^PfX{x g )f, (20) 

/ 

where the sum is taken over all Fock states / in the r\ c meson, pf is the probability to find the Fock state /, (2:3) / is 
the fraction of momentum carried by gluon in the Fock state /. Assuming that the dominant contribution to formula 
(|20| is given by the quark-antiquark-gluon Fock state one gets (x g ) n ~ 3 • 10~ 4 . 

Further let us pay attention to formula fj 15[) . This formula allows us to calculate the mean energy of the gluon in 
the rj c meson 

(77 e |^+a J ((iV )gg J )x|0) , f/ > / nfi(: xni2^v ai\ 
E g = — 7—s ' — = M v {x 3 ) = (0.66 ± 0.18) GeV. (21) 

(Vc\^ + cr ] {gH j )x\0} 

One can also estimate the gluon energy from formula (TTBl) 



(v^a-wl - t) 9 H 3 )m = [h^ = 0M GeV . (22) 



\ M^+a^gH^xlO) 



This value is in a good agreement with (|2"Tj) . 

In the last excise the relative momentum p r of the color-octet quark-antiquark pair will be determined. To do this 
let us apply formula (fTT| 



1 



_2 ^2 



(r 1c \TP+a j ((-V k gH j -gH j V k + 2V k gH j V k ))x\a) _ 1 l((x 1 -x 2 f) 



= oW 3 >\ / VJ M v = 0.49 GeV, (23) 



IPr| ~ 2(1-^3))^ {r ]c \^+a 3 {gH ] )x\0) 2\) (1 - (x 3 )) 

It is important to note that the relative velocity of the color-octet quark-antiquark pair {v 2 ~ 3((xi— a:2) 2 )/(l— (^3)) 2 ~ 
0.11) is two times smaller than that for the color-singlet quark-antiquark state(i) 2 ~ 0.21 Q). What is expected result 
since the gluon carries some fraction of momentum and energy of the r\ c meson. 



V. THE MODEL OF THE DISTRIBUTION AMPLITUDE V{xi,x 2 ,x 3 ). 

To build the model of the DA V(xi. X2, X3) one can recall that the models of the leading twist charmonia distribution 
amplitudes proposed in papers 0, 0, fll| can be written in the following form 

/ M 2 -\ 

^-^•expf-^J, (24) 



5 



where 4> as is the asymptotic form of the distribution amplitude <j>, M 2 S is the invariant mass of the cc system expressed 
in terms of kinematic variables, M 2 is the characteristic invariant mass of the cc system. 

In the case of the quark- antiquark-gluon Fock state the r\ c meson can be divided into two subsystems: the quark- 
antiquark pair in the color-octet state and the color-octet pair with the gluon. If one introduces the following variable 

xi = xy, x 2 = (l-x)y, x 3 = l-y, (25) 

the invariant masses of the corresponding subsystems can be written in the form 



Ml = ml— i_ = m 2iL^L, M { 2 cE}g =mg 1 =m giL-g!). (26) 
x(l - x) x\x 2 1 19 xy(l-x) x\x 2 



To build the model of the DA V(xi, x 2 , X3) one should introduce the exponent factors for the both subsystems. Thus 
one has 

V(x 1: x 2 , x 3 ) = c(p 1 ,p 2 )xix 2 xl ■ exp f-ft. (1 ~ X ^ ) ■expf- / g 2 (1 ~ a:3) Y (27) 

V XiX 2 J V X 1 X 2 J 

where the constant c(/3i, fi 2 ) can be determined from normalization condition @. The values of the parameters /3i, j3 2 
can be fixed from the requirement that the values of the moments (X3), ((x\ — x 2 ) 2 ) for model (I2T[) must coincide with 
the values (JT5J, CE3]). Thus one gets ft = 1.25±g'||, /3 2 = 1.15 ± 0.25. Note that model (J27J) is defined at the scale 
(X ~ m c . 

One can assume that it is sufficient to introduce only one exponent with the total invariant mass of the quark- 
antiquark-gluon system in formula (|27)l . However, the calculation shows that in this case it is not possible to reproduce 
values (fT2|) . (|13|) of the moments (2:3), ((xi — x 2 ) 2 ) simultaneously. Physically, this means that although energy scales 
of the two subsystems introduced above are of the same order they are not equal to each other. 



VI. CONCLUSION 



In this paper the properties of the quark-antiquark-gluon Fock state in the r\ c meson are studied. This properties 
can be parameterized by the moments of the twist-3 distribution amplitude ([1]) and the constant f^. To calculate 
the first, the second order moments and the constant QCD sum rules was applied. The result of this calculation 
allowed one to build the model of the DA which can be used in the study of different production processes. 

To make the results of the calculation more transparent, QCD operators which determine the moments of the 
distribution amplitude are expanded in relative velocity of the quark-antiquark pair in the rj c meson. In particular, 
it was shown that the quark-antiquark pair is in the S'-wave spin-triplet color-octet state. The probability amplitude 
to find the quark-antiquark pair in the color-octet state is rather small ~ 0.04. The energy of the gluon in the quark- 
antiquark-gluon Fock state is ~ 600 MeV. The relative velocity of the color-octet quark-antiquark pair is ~ 0.11 what 
is much smaller than that in the color-singlet quark-antiquark pair. 
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